Dynamical Subgrid-scale Parameterizations for Quasigeostrophic Flows using Direct Numerical Simulations

In this thesis, parameterizations of non-linear interactions in quasigeostrophic (QG) flows for severely truncated models (STM) and Large Eddy Simulations (LES) are studied. Firstly, using Direct Numerical Simulations (DNS), atmospheric barotropic flows over topography are examined, and it is established that such flows exhibit multiple equilibrium states for a wide range of parameters. A STM is then constructed, consisting of the large scale zonal flow and a topographic mode. It is shown that, qualitatively, this system behaves similarly to the DNS as far as the interaction between the zonal flow and topography is concerned, and, in particular, exhibits multiple equilibrium states. By fitting the analytical form of the topographic stationary wave amplitude, obtained from the STM, to the results obtained from DNS, renormalized dissipation and rotation parameters are obtained. The usage of renormalized parameters in the STM results in better quantitative agreement with the DNS. In the second type of problem, subgrid-scale parameterizations in LES are investigated with both atmospheric and oceanic parameters. This is in the context of two-level QG flows on the sphere, mostly, but not exclusively, employing a spherical harmonic triangular truncation at wavenumber 63 (T63) or higher. The methodology that is used is spectral, and is motivated by the stochastic representation of statistical closure theory, with the ‘damping’ and forcing covariance, representing backscatter, determined from the statistics of DNS. The damping and forcing covariance are formulated as 2 × 2 matrices for each wavenumber. As well as the transient subgrid tendency, the mean subgrid tendency is needed in the LES when the energy injection region is unresolved; this is also calculated from the statistics of the DNS. For comparison, a deterministic parameterization scheme consisting of 2× 2 ‘damping’ parameters, which are calculated from the statistics of DNS, has been constructed. The main difference between atmospheric and oceanic flows, in this thesis, is that the atmospheric LES completely resolves the deformation scale, the energy and enstrophy injection region, and the truncation scale is spectrally distant from it, being well in the enstrophy cascade inertial range. In oceanic flows, however, the truncation scale is in the vicinity of the injection scale, at least for the parameters chosen, and is therefore not in an inertial range. A lower resolution oceanic LES at T15 is also examined, in which case the injection region is not resolved at all. For atmospheric flows, it is found that, at T63, the matrix parameters are practically diagonal so that stratified atmospheric flows at these resolutions may be treated as uncoupled layers as far as subgrid-scale parameterizations are concerned. It is also found that the damping parameters are relatively independent of the (vertical) level, but the backscatter parameters are proportional to the subgrid flux in a given level. The stochastic and deterministic parameterization schemes give comparably good results relative to the DNS. For oceanic flows, it is found that the full matrix structure of the parameters must be used. Furthermore, it is found that there is a strong injection of barotropic energy from the subgrid scales, due to the unresolved, or partially resolved, baroclinic instability injection scales. It is found that the deterministic parameterization is too numerically unstable to be of use in the LES, and instead the stochastic parameterization must be used